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ABSTRACT 

We present a measurement of the distant Type la supernova rate derived 
from the first two years of the Canada - France - Hawaii Telescope Supernova 
Legacy Survey. We observed four one-square degree fields with a typical temporal 
frequency of {At) ~ 4 observer-frame days over time spans of from 158 to 211 days 
per season for each field, with breaks during full moon. We used 8-10 meter-class 
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telescopes for spectroscopic followup to confirm our candidates and determine 
their redshifts. Our starting sample consists of 73 spectroscopically verified Type 
la supernovae in the redshift range 0.2 < z < 0.6. We derive a volumetric SN la 
rate of r v ((z) = 0.47) = 0.42+°;^ (systematic) ±0.06 (statistical) xl0~ 4 yr" 1 
Mpc 3 , assuming h = 0.7, Q m = 0.3 and a flat cosmology. Using recently published 
galaxy luminosity functions derived in our redshift range, we derive a SN la rate 
per unit luminosity of r L ((z) = 0.47) = O.^^;^ (systematic) t^oii (statistical) 
SNu. Using our rate alone, we place an upper limit on the component of SN la 
production that tracks the cosmic star formation history of 1 SN la per 10 3 M Q of 
stars formed. Our rate and other rates from surveys using spectroscopic sample 
confirmation display only a modest evolution out to z — 0.55. 

Subject headings: galaxies: evolution - galaxies: high redshift - supernovae: 
general 



1. Introduction 



Type la supernovae (SNe la) have achieved enormous importance as cosmological dis- 
tance indicators and have provided the first direct evidence for the dark energy th at is driving 
the Universe's accelerated expansion (IRiess et al.lll998l : iPerlmutter et al.lll999l ). In spite of 
this importance, the physics that makes them such useful cosmological probes is only partly 
constrained. White dwarf physics is the best candidate for producing a standard explosion 



1 Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope (CFHT) which is operated by the National Re- 
search Council (NRC) of Canada, the Institut National des Sciences de 1'Univers of the Centre National 
de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part 
on data products produced at the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii 
Telescope Legacy Survey, a collaborative project of NRC and CNRS. Based on observations obtained at 
the European Southern Observatory using the Very Large Telescope on the Cerro Paranal (ESO Large Pro- 
gramme 171.A-0486). Based on observations (programs GN-2004A-Q-19, GS-2004A-Q-11, GN-2003B-Q-9, 
and GS-2003B-Q-8) obtained at the Gemini Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini part- 
nership: the National Science Foundation (United States), the Particle Physics and Astronomy Research 
Council (United Kingdom), the National Research Council (Canada), CONIC YT (Chile), the Australian 
Research Council (Australia), CNPq (Brazil) and CONICET (Argentina). Based on observations obtained 
at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute 
of Technology, the University of California and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support of the W.M. Keck Foundation. 
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due to the well understood Chandrasekhar mass limit (1 Chandr asekharl 1 1 93 ll ) . However, any 
plausible SN la scenario requires a companion to donate mass and push a sub-Chandras ekhar 
C-0 white dwarf towards this limit producing some form of explosion (for a review, see iLivio 
200ll ). The range of possible companion scenarios needed to accomplish this are currently 
divided into two broad categories: the single degenerate scenario, where the companion 
i s a subgiant or giant star that is donating matter through winds or Roche lobe overflow 
rtwhelan fc l"ben]| 19731 : lNomotd[l~982l : banal et al.lll996l : Irian fc Podsiadlowski|[20ol , and the 
double degenerate scenario involving the coalescence of tw o white dwarf stars after losing 
orbital angular momentum through gravitational radiation (jWebbinklll984l ; llben &: Tutukov 
19841 : iTornambe fc Matteuccilll986l : iNapiwotzki et al.ll2004h . 



Population synthesis models for these sce narios predict different SN la production 



timescales relative to input star formation (e.g. iGreggiol 120051 ). By comparing the global 
rate of occurrence of SN la at different redshifts to measurements of the global cosmic star 
formation history (SFH), the 'delay function', parameterized by its characteristic timescale, 
r, can be derived, which in turn constrains the companion scenarios. This comparison re- 
quires the calculation of a volumetric SN la rate and measuring the evolution of this rate 
with redshift. 



Early SN surveys were host targeted (e.g., IZwickyill938l ). and produced rates per unit 
blue luminosity that required conversion to volumetric rates through galaxy luminosity func- 
tions. These surveys suffer from large systematic uncertainties because of the natural ten- 
dency to sample the brighter end of the host luminosity function. With the advent of wide- 
field imagers on moderately large telescopes, recent surveys have been able to target specific 
volumes of space and directly calculate the volumetric rate. Examples of volumetric SN la 
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We plot the ra tes from these surveys as a function of redshift in Figure [U, along with a 
recent SFH fit from lHopkins fc Beacoml (120061 ). renormalized by a factor of 10 3 . This allows 
us to compare the SFH to the observed trend in the SN la rate. This trend, compared 
with the SF H curve, shows some cu rious properties. The large gradient just beyond z = 0.5 
observed in iBarris fc Tonryl (120061 ) has no analog in the SFH curve, a nd neither does the 
apparent down-turn beyond z = 1.2 observed by iDahlen et al.l (120041 ) . Fits of the delay 
function to various subsets of these data have produced no consensu s on r, or the form o f 
the delay. Reported values f or r range from as s hort as r < 1 Gyr (IBarris fc Tonryl 120061 ) 
to as long as r = 2 — 4 Gyr (jStrolger et al.l 120041 ) . This lack of consensus and the peculiar 
features in Figure [1] argue that systematics are playing a role in the observed SN la rates, 



-4- 



especially at higher redshifts. It is vital to investigate the sources of systematic error in 
deriving SN la rates and to compare the cosmic SFH with rates that have well characterized 
systematic errors. 

In this paper we take advantage of the high-qualit y spectroscopy an d well-defined sur- 



vey properties of the Supernova Legacy Survey (SNLS. lAstier et al.ll2006l ) to produce a rate 
that minimizes systematics, which we then compare with cosmic SFH. To minimize contam- 
ination, we use only spectroscopically verified SNe la in our sample. We examine sources 
of systematic error in detail and, using Monte Carlo efficiency experiments, place limits on 
them. In particular, we improve upon previous surveys in the treatment of host extinction 



by using the recent dust models of iRiello fc Patatl (120051 ). We also investigate the possibility 



that SNe la are being missed in the cores of galaxies with fake SN experiments using real 
SNLS images. These experiments allow us to place limits on our own errors and assess the 
impact of various sources of systematic error on SN surveys in general. 

In order to avoid large and uncertain completeness corrections, we employ a simpli- 
fication in our rate determination that fully exploits the data set currently available: we 
restrict our sample to the redshift range 0.2 < z < 0.6. This ensures that the majority 
of SNe la peak above our nominal detection limits and thus provides a high completeness. 
This simplification allows us to extract a well-defined sample of spectroscopically verified 
SNe la from the survey and accurately simulate the SNLS survey efficiency, thus producing 
the most accurate SN la rate at any redshift. Our rate alone is sufficient to constrain some 
of the SFH delay function models by placing limits on their parameters. 

The structure of the paper is as follows. In §|2]we describe the survey properties relevant 
to rate calculation. In §|3] we develop objective selection criteria, derive our SN la sample 
and analyze this sample to determine our spectroscopic completeness. In §3]we describe our 
method for calculating the survey efficiency and present the results of these calculations. In 
§5] we present the derived SN la rates per unit volume and per unit luminosity, an analysis 
of systematic errors, and a comparison of our rates with rates in the literature. In §6] we 
compare our volumetric rate and a selection of rates from the literature with two recent 
models connecting SFH with SN la production. 

For ease of comparison with other rates studies in the literature, we assume a flat 
cosmology throughout with H = 70 km s _1 Mpc -1 , = 0.7, and — 0.3. 
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2. The Supernova Legacy Survey 



The SNLS is a second-generation SN la survey spanning five years, instigated with the 
purpose of measuring the accelerated expansion of the universe and constraining th e average 
pressure-density ratio of the universe, (w), to better than ±0.05 (lAstier et al.ll2006l ). In order 
to achieve this goal, the SNe la plotted on our Hubble diagram must have well-sampled light 
curves (LCs) and spectral followup observations that provide a ccurate redshifts an d solid 
identifications. The LC sampling is achieved using MegaCam (IBoulade et al.ll2003l ). a 36 
CCD mosaic one-square degree imager, in queued service observing mode on the 3.6 meter 
Canada- France-Hawaii Telescope (CFHT). This combination images four one-square degree 
fields (Dl-4, evenly spaced in right ascension, see Table [T]) in four filters (g'r'i'z 1 ) with an 
observer-frame cadence of At ~ 4 days (rest-frame cadence for a typical SN of At ~ 3 days) 
and with a typical limiting magnitude of 24.5 in i'. The queued service mode provides robust 
protection against bad weather, as any night lost is re-queued for the following night. 

This observing strategy provides dense LC coverage for SNe la out to z ~ 1 and is ideal 
for measuring the rate of occurrence of distant SNe la. It also produces high quality SN la 
candidates identified early enough so that spectroscopic follow up observations can be sched- 
uled near t he candidate's max imum light (jSullivan et al.ll2006l ). This strategy has been very 
successful (IHowell et al.ll2005l ). and the SNLS has been fortunate to have consistent access to 
8-10 meter class telescopes (Gemini, Keck, VLT) for spectroscopic followup. This is critical 
for providing a high spectroscopic completeness and the solid spectroscopic t ype confirma- 



tion required to remo ve contaminating non-SN la objects from our sample (IHowell et al. 
2005l : lBasa et al.ll2006h . 



2.1. The Detection Pipeline 

The imaging data are analyzed by two independent search pipelines in Canada^ and 
Franc^E For the rate calculation in this paper, we use the properties of the Canadian 
pipeline. 

The Canadian SNLS real-time pipeline uses the i' filter images for detection of SN 
candidates and images in all filters for object classification. Each epoch consists of five to 
ten exposures which undergo a preliminary (real-time) reduction which includes a photo- 
metric and astrometric calibration before being combined. A reference image for each field 



2 see |http: / /legacy.astro.utoronto.ca/l 



see http: / /makiki.cfht .hawaii.edu:872 /sne / 
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is constructed from previously acquired, hand picked, high-quality images. The detection 
pipeline then seeing-matches the reference image to the (usually lower image quality) new 
epoch image (IPritchet et al.ll2006l ). The seeing-matched reference image is then subtracted 
from the new epoch and the resulting difference image is analyzed to detect variable objects 
which appear as residual (positive) point sources. A final list of candidate variable objects 
is produced from this difference image in two stages: first, a preliminary candidate list is 
generated using an automated detection routine and then, a final candidate list is culled by 
human review of the preliminary list. This visual inspection is conducted by one of us (D.B.) 
and is essential for weeding out the large quantity of non-variable objects (image defects, 
and PSF matching errors) that remain after the automated detection stage. 

At this stage, all candidate variables are given a preliminary classification and any object 
that may possibly be a SN (of any type) has 'SN' in its classification. The new measurements 
of variable candidates are entered into our object database and compared with previously 
discovered variable objects. This comparison weeds out previously discovered non-SN vari- 
ables such as AGN and variable stars from the SN candidate list. All measurements of the 
current SN candidates, including recent non-detections, are then evaluated for spectroscopic 
followup using photometric selection criteria. 



The details of the photometric selection process for the SNLS are presented in lSullivan et al. 



(120061 ). In brief, all photometric observations of the early part of the LC of a SN candidate 
are fit to template SN la LCs using a \ 2 minimization in a multi-parameter space that in- 
cludes redshift, stretch, time of maximum light, host extinction, and peak dispersion. The 
template LCs are generated f rom an updated version of the SN la spectral templates pre- 
sented in iNugent et al.l (120021 ) . These spectral templates are multiplied by the MegaCam 



filter response functions and integrated, thus accounting for k-corrections ([Sullivan et al. 



20061 ). The results of this fit are used to measure a photometric redshift, z PHO t, for the 
candidate and to make a more accurate classification. If there is any doubt about the nature 
of the object, the 'SN' classification is retained in the database. 

All SN candidates in the database are available for the observers doing spectroscopic 
followup. The quality of the candidate, deduced from the template fit and an assessment 
of usefulness for cosmology, is used to prioritize the candidates for spectroscopic observa- 
tion. Once these observations are taken, they are reduced and compared to SN la spectral 
templates (IHowell et al.ll2005l ; iBasa et al.ll2006l ) to calculate a spectroscopic redshift, zspec- 
The final typing assessment uses all available information, both photometric and spectro- 
scopic. The photometry provides early epoch colors, which can help identify CC SNe. It 
also provides an accurate phase for the spectroscopic observation which is also important 
in discriminating SNe la from CC SNe. The galaxy-subtracted candidate spectrum is then 
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checked for the presence of spectral features peculiar to SNe la. We then assig n a likelihood 



statistic for the candidate's membership in the SN la type (IHowell et al.ll2005l ) 



3. Selection Criteria 

Selection criteria are used to provide consistency between the observed sample, the sur- 
vey efficiency calculation, and the completeness calculation and thus produce an accurate 
rate. In practice, they serve to objectify the survey goals and properties (which unavoidably 
include the human element) such that efficiency simulations are accurate and tractable. The 
criteria we developed consist of the minimum required photometric observations, expressed 
in terms of rest-frame epoch and filter, that guarantee that any real SN la acquires spec- 
troscopic followup. They were derived by examining the photometric observations of all 
our spectroscopically confirmed SNe la in the redshift range 0.2 < z < 0.6. To account 
for any real SN la that meet these criteria but, for one reason or another, did not acquire 
spectroscopic followup, we also apply these criteria to our entire 'SN' candidate list in a 
completeness study (see below). 

Since the primary goal of the SNLS is cosmology, when selecting SN candidates for 
spectroscopic followup we attempt to eliminate objects, even SNe la, that offer no information 
for cosmological fitting. Examples of these include SNe for which no maximum brightness 
can be determined, or for which no stretch or no color information can be measured. Thus, 
the objective criteria that define our sample and survey efficiencies are expressed by requiring 
each confirmed SN la to have the following observations: 

1. one i' detection at S/N > 10.0 between restframe day -15.0 and day -1.5 

2. two i' observations between restframe day -15.0 and day -1.5 

3. one r' observation between restframe day -15.0 and day -1.5 

4. one g' observation between restframe day -15.0 and day +5.0 

5. one i' or r' observation between restframe day +11.5 and +30.0 

Criterion 1 and 2 implement our need to detect candidate SNe la early enough to sched- 
ule spectroscopic observations near maximum brightness. Criterion 2 is required to judge if 
the LC is rising or declining. Criteria 3 and 4 are required because an early color is important 
for photometrically classifying the SN type. Criterion 5 implements the requirement that 
stretch information be available for any cosmologically useful SN la. We only require the 
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detection in the pre-max %' because during the early part of the light curve SNe la are dis- 
tinguished from other SN types by having redder colors. Thus, if we have an early detection 
of a candidate in i', but can only place a limit on the object in r' or g', then it would have a 
reasonably high probability of being a SN la and is likely to be spectroscopically followed up. 
This also means that highly reddened SNe are not selected against. For this redshift range, 
we need not be concerned with criteria based on the z' filter. Criterion 5 would not logically 
enter the selection process as a detection, since these observations could not be taken before 
the decision to followup is made. It is included solely to remove objects that are discovered 
close to the end of an observing season, when there is no hope of obtaining the observations 
needed to derive a stretch value. 

It is important to point out that these criteria are independent of the LC fitting that 
is normally done in candidate selection, i.e., there are no criteria involving the SN la fit \ 2 ■ 
This is because we defined these criteria with spectroscopically confirmed SNe la. The fitting 
is required to derive the type and the redshift of candidate SNe. In our sample selection 
both of these quantities are given by the spectroscopy. In the efficiency simulations, we are 
only interested in our detection efficiency for SNe la, so the type is defined a priori, and the 
redshifts are given by the Monte Carlo simulation (see below). The LC fitting does enter 
into the completeness study, since we are then interested in objects without spectroscopy. 
We describe the LC fitting criteria used to derive an accurate completeness, given the above 
selection criteria, in §3.21 



3.1. The Observed Sample 

In order to define an observed sample consistent with these survey selection criteria, we 
must eliminate spectroscopically confirmed SNe la in the initial list that do not meet these 
criteria. These 'special-case' SNe la acquired spectroscopic followup for two reasons. First, 
during some of our initial runs we attempted to spectrally follow up nearly every suspected 
SN to help refine our photometric selection criteria. Second, occasionally bad weather can 
prematurely end a field's observing season before all the good declining candidates have the 
required observations to determine their stretch values. 

We derived our starting SN la sample from all spectroscopically confirmed SNe la with a 
spectroscopic redshift, z S p EC , in the range 0.2 < zspec < 0.6, discovered in the first two full 
seasons of each deep field. The starting and ending dates and resulting time span in days is 
listed for each season of each field in Table [TJ Figure [3] illustrates the field observing seasons 
for the sample by plotting the Julian Day of the epochs versus their calculated limiting i' 
magnitude (see §4.ip . 
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Table [2] individually lists the 73 spectroscopically confirmed SNe la from the SNLS that 
comprise our starting sample. Column 1 gives the SNLS designation for the SN, columns 
2 and 3 give the J2000.0 coordinates, column 4 gives the redshift, column 5 gives the MJD 
of discovery, column 6 indicates if the object was culled from the initial list by enumerating 
which of the criteria from $3] it failed, and column 7 lists the references in which further 
information about the object is published. The table is ordered by field and by time within 
each field with breaks between the two seasons of the given field. 

Figure [2] plots the nightly averaged photometry for each of the 73 spectroscopically 
confirmed SNe la from Table [2] with la err or bars on a norma lized AB magnitude scale. 
The best-fit SN la template is overplotted (jSullivan et al.ll2006l ). The magnitude scale is 
normalized such that the brightest tick mark is always 20 magnitudes. This procedure 
preserves the relative magnitude difference between filters for a given SN. The day scale on 
the bottom of each plot is the observed day relative to maximum light. The day scale on the 
top is the restframe day relative to maximum light. The designation from Table [2] (minus 
the SNLS prefix) is given in the upper left corner and the spectroscopic redshift in the upper 
right corner of each panel. If the object was culled from the initial list, this is indicated 
under the designation with the word 'Rej' (e.g. SNLS-03Dldj was culled). 

Almost every object from season one of each field is included in the cosmology fit of 
Astier et al.l (120061 ). The three exceptions are SNLS-03Dlar, which had insufficient obser- 
vations at the time, SNLS-03D4cj, which was a SN 1991T-like SN la, and SNLS-03D4au, 
which was under-luminous, most likely due to extinction. We point out that even though 
these objects were excluded from the cosmology fit, their identity as SNe la has never been 
in doubt. The obj ects that were observed with Gemini have their spectra published in 
Howell et al.l (2005 ). The o bjects observed with the VLT will have their spectra published 
shortly in iBasa et al.l (120061 ). The remaining 10 objects from the first seasons were observed 
with Keck. 

The sample is summarized in Tabled which lists, for each season of each field, the total 
number of spectroscopically confirmed SNe la and the number after culling the starting list 
using our objective selection criteria. 



3.2. Spectroscopic Completeness 

We now calculate the number of objects that passed our selection criteria but, for one 
reason or another, were not spectrally followed up. This calculation is aided by the high 
detection completeness of the survey below z = 0.6 (see §4.2.31) . and the classification scheme 
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we use, where any object remotely consistent with a SN LC, after checking for long-term 
variability, retains the 'SN' in the classification. We are also able to use a final version 
of the photometry, generated for all objects in our database from images that have been 
de-trended with the final calibration images for each observing run. This final photometry, 
which now covers all phases of the candidate LCs, is fit with SN la templates as described 
above to produce a more accurate zphot and a XsNia f° r the \ 2 °f the SN la template fit 
to the photometry in all filters. We examined all objects with a final photometry Zphot 
in the range 0.2 < zphot < 0.6 discovered within the time spans in Table [T] with the 
following classifications: 'SN', 'SN?', 'SNF, 'SNIP, 'SNII?', 'SN/AGN', and 'SN/var?'. We 
measured the offset and uncertainty in our z PHO t fitting technique by comparing zphot 
with zspec and found a mean offset of Az < 10~ 3 and an RMS scatter of a z = 0.08. We, 
therefore, assume that the remaining error in Zphot from the the final photometry is small 
and random such that as many candidates are scattered out of our redshift range of interest 
as are scattered in. 

Of 180 objects from the sample time ranges with 'SN' in their type, 50 do not have the 
required observations from our object selection criteria listed above and so, even if they were 
SNe la, would not be included in our culled sample. Of the remaining 130 objects, 64 are 
rejected because their fit to the templates has a XsNia > 10-0 anc ^ so ver y unlikely to be 
SNe la rtSullivan et aDbood ). We then apply an upper limit stretch cut, requiring s < 1.35, 



to the remaining 66 objects. Objects with s > 1.35 are also not SNe la (see Astier et al. 
2005, Figure 7 and Sullivan et al. 2006b, Figure 3). These objects are probably SNe IIP 
which have a long plateau in their LCs and hence produce anomalously high s values when 
fit with a SN la template. The s < 1.35 cut removes 33 objects. We then make a cut by 
examining the early colors and remove those that have large residuals in this part of the 
LC as a result of being too blue (one signature of a core-collapse SN). Of the remaining 33 
objects, 14 of these are rejected as too blue in the early colors, even though the overall XsNia 
is less than 10.0. We are left with 19 unconfirmed SN la candidates that have a reasonable 
probability of being missed, real SNe la. 

Table 0] lists the 19 unconfirmed SN la candidates, their coordinates, their zphot, 
discovery date, initial type, XsNiai anc ^ their status. We group them into those with XsNia < 
5.0 and those with XsNia > 5-0 and consider those in the first group to be probable SNe la, 
and those in the second group to be possible SNe la. We point out that the intrinsic variation 
in SN la LCs rarely allow templat e fits with XsNia < 2, and that typical fits have XsNia i n 
the ranee 2-3 JSullivan et alJbood ). We take the conservative approach that, aside from the 



division at X%Nia = 5-0, we must consider each candidate in each group as equal. This then 
determines the range of completeness we consider in calculating our systematic errors (see 
below). Our most likely completeness is defined by assuming that each 'Probable SN la' in 
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this list is, in fact, a real SN la and each 'Possible SN la' is not. The minimum completeness 
is defined by the scenario that all 19 are real SN la and the maximum completeness is 
defined by the scenario that none of the 19 are real, which amounts to 100% completeness. 
We tabulate the confirmed, probable and possible SNe la and the minimum and most likely 
completeness for each field and the ensemble in Table [5j We will use this table when we 
compute our systematic errors in §5.3.11 

Figure H] plots the nightly average photometry for the 19 unconfirmed SN la candidates 
from Table H] using the same normalized AB magnitude scale and day axes as in Figure [2J 
The photometric redshift is indicated in the upper right corner of each panel. The XsNia 
value from Table H] for each SN is indicated under its designation on each panel. 



4. Survey Efficiency 



Since Fritz Zwicky's pioneering eff orts to estima te supernova rates from photographic 
surveys using the control-time method (jZwickvlll938l ). there ha ye been significant improve- 
ments in calculating a give n survey's efficien cy (for a review see IWood-Vaseyll2005l . Chapter 
6). As a recent example, iPain et al. Jl996h used SN la template LCs to place simulated 
SNe in CCD survey images to generate a Monte Carlo simulation that produced a much 
more accurate efficiency for their survey. Most recent surveys using CCDs have performed 
some variation of this method to calculate their efficiencies and from them derive their rates 



flHardin et al.ll2000l : IPain et all 12002 ; Madgwick et al.ll2003l : Blanc et alJl2004f ). 



In our particular variation on this method, we do not place artificial SNe on every image 
of our survey. Instead, we characterize how our frame limits vary with relevant parameters 
(such a seeing) using a subset of real survey images. We then use this characterization 
to observe a Monte Carlo simulation which uses the updated SN la spectral templates of 
Nugent et al.l (120021 ) and our survey filter response functions to generate the LCs from a 
large population of realistic SNe la. Thus, to calculate an appropriate survey efficiency, we 
need to implement the objective selection criteria defined above in a Monte Carlo efficiency 
experiment that simulates the observation of the SN la LCs by the SNLS. Criteria 2-5 (see 
£J3]) can be implemented simply by inputting the date and filter of each image in the survey 
sample time ranges, and seeing if we have the required observations for each simulated 
candidate. Criterion 1 specifies a detection in the i' filter, which requires that we calculate 
the SN visibility at each %' epoch in the survey sample time ranges. 
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4.1. %' SN Visibility 

The photometric depth reached by a given %' observation depends on the exposure 
time (E e ), image quality (IQ e ), airmass {X e ), transparency (T e ), and the noise in the sky 
background {S e ). Some of these data are trivially available from each image header. The 
transparency and the sky background must be derived from the images themselves. 

Our final photometry pipeline includes a photometric calibration process that calculates 
a flux scaling parameter, F e , for each image. We calculate it by comparing a large number of 
isolated sources in the object image with the same objects in a (photometric) reference image. 
The resulting F e values are applied to each object image to ensure that the flux measured 
for a non-variable object is the same in each epoch. Thus, F e accounts for variations in both 
T e and X e . An image with lower transparency and/or higher airmass will have a larger F e . 
During this process the standard deviation per pixel in the sky is also calculated, allowing us 
to account for variations in S e . The total number of usable CCD chips, out of the nominal 
36, is also tabulated (see §4.1.21) . 

Another factor that determines a spatially localized frame limit is the galaxy host back- 
ground light against which the SN must be discerned (Hi^ ga i). This depends on the brightness 
and light profile of the host and the brightness and position of the SN within the host. This 
dependence is mitigated somewhat by the subtraction method used in our detection pipeline 
(see §2.ip . but must still be measured. 

We designed a controlled experiment to explore the effects of IQ e , and FLi^ ga i on SN 
visibility. This experiment places many artificial SNe of varying brightness and host galaxy 
position (yielding a range of Hi^ ga i) in real SNLS detection pipeline images of varying IQ e . 
We chose a range of IQ e from IQ e = 0'.'60, close to the median for the survey, to IQ e = 1'.'06, 
near the limit of acceptability. We used epochs with the canonical exposure time of 3641s 
and required that the images were taken under photometric conditions. 



Prior to the addition of fake SNe, each image was analyzed with SExtractor (IBertin fe Arnouts 



19961 ) to produce a list of potential galaxy hosts over the entire image. For a given fake SN, 
the host was chosen from this list using a brightness weighted probability, such that brighter 
galaxies are more likely to be the host than fainter galaxies. The location within the host 
for the fake SN was also chosen with a brightness weighted probability, such that more SNe 
are produced where the galaxy has more light (i.e. toward the center). Once the location 
within the pipeline image is decided, a nearby isolated, high S/N star was scaled to have a 
magnitude in the range 21.0 < i' < 27.0 and added at the chosen position. 



There was no correlation of the fake SN magnitude with the host magnitude, therefore, 
our simulations were relevant for SNe at all phases of their LC. This spatial distribution and 
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magnitude range allows us to quantify any systematic loss of SN visibility near the cores of 
galaxies in our recovery experiments (see below). Once a set of these images was produced 
it was pu t through the same detection pipeline used by the Canad ian SNLS for detecting 



real SNe flPerrett et al.ll2006l ; ISullivan et al.ll2006l : lAstier et al.l 120061 ) 



Figure [5] shows the raw recovery percentage of ~ 2000 fake SNe after the human review 
process for two IQ e values: 0'.'69 and 1'.'06. The 50% recovery limits are indicated and 
are the most useful for rate calculation since the visible SNe missed below these limits are 
gained back by including the invisible SNe above the limits (see Figure Eh)- The loss in 
visibility going from automatic detection to human review amounts to a brightening of the 
visibility limits of only 0.1 magnitudes at the small IQ e value. Figure [5b shows no trend 
with host offset and Figure EH shows that the cutoff due to background brightness is 20 
mag arcsecond -2 . A notable feature of Figure Eh is the maximum recovery percentage of 
95% for the IQ e = 0'.'69 image. We examined the spatial distribution of the fake SNe from 
this image that were missed above i' = 23 to try to understand the source of this limit on 
the recovery. We saw no correlation with galaxy host offset, proximity to bright stars, or 
placement on masked or edge regions. This feature appears to be purely statistical in origin 
and we account for it when observing the Monte Carlo simulations (see §4.2p . 

Figure [6] shows the 50% recovery limits derived from the fake SN experiments using nine 
i' images having a range of IQ e - These limits have been corrected for sky noise, transparency, 
and exposure time differences. We plot the histogram of all the IQ e values for all i' images 
relevant to this study as a dashed line. All points are derived at the automated detection 
stage unless otherwise indicated. The corrected results of the human review experiment from 
above are plotted as asterisks. We fit the human review limits with a linear fit (shown as the 
solid line) and this fit represents an upper bound on the limits of the images we sampled. 
A constant frame limit of i' = 24.5 is shown as the solid horizontal line and is a reasonable 
lower bound. The two solid lines encompass all points in Figure El We will use the human 
review limit fit as our best estimate of the frame limit versus IQ e function, with the constant 
limit as an estimate for the systematic error in our rates due to the %' frame limits (see §5.3p . 



4-. 1.1. %' SN Visibility Equation 

Our fake SN experiments have provided a way to calculate the visibility limit in mag- 
nitudes, L e , for any i' epoch in our survey sample time span using the following formula: 

L e = L .5 - a(IQ e ~ 0.5) + 2.5 log{E e /E ref ) - 2.5 log(F e ) - 2.5 \og(S e / S ref ), (1) 
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where L0.5 is reference visibility limit for an epoch with IQ e = 0"5, F e = 1.0, exposure time 
of E re f seconds and sky noise of S re f counts, a is the proportionality factor between IQ e 
and the visibility limit, E e is the exposure time of the epoch, F e is the flux scale factor, and 
S e is the sky noise in counts of the epoch (see Figure [3]). This formula assumes a linear 
relationship between IQ e and L e , which appears to be a reasonable approximation over the 
range of IQ e used to discover SNe (see Figure [6]). 

Table lists the parameters calculated using the 50% recovery fraction visibility limits 
determined from the human review recovery experiment (see Figure and [H]) • Columns 1 
and 2 of Table [6] lists the IQ e for the pair of good and bad IQ images used in the human 
review experiment, column 3 lists the reference exposure time, column 4 lists the reference 
sky noise in counts, column 5 lists the visibility limit at IQ = 0'.'5, and column 6 lists the 
proportionality constant between IQ e and L e . For the reference sky noise, S re f, we used the 
value from the good IQ image and adjusted the limit from the poor IQ image to correspond 
to an image with the same sky noise as the good IQ image. 

4-1-2. Temporary CCD Losses 

Another factor affecting the visibility of SNe in the SNLS must be accounted for. Oc- 
casionally, a very small subset of the 36 MegaCam CCDs will malfunction for a short time, 
usually because of a failure in the readout electronics. An even rarer occurrence is the 
appearance of a condensate of water on the surface of one of the correctors that covers a 
localized area of the field of view rendering that part of the detector temporarily useless for 
the detection of SNe. When we calculate the fluxscale factors mentioned above, the num- 
ber of usable CCDs is also recorded. This number is used to account for these localized, 
temporary losses of SN visibility (see below). 

4.2. Monte Carlo Simulation 

The Monte Carlo technique allows us to determine our survey efficiency to a much higher 
precision than permitted by the small number of observed events. Using observed SN la LC 
properties and random number generators, we simulate a large (N = 10 6 ) population of SN la 
events in the sample volume occurring over a two year period centered on the observed 
seasons for the field. This large number is sufficient to drive the Poisson errors down to 
vN /N = 0.1%. This population is then observed by using real SNLS epoch properties and 
equation [U combined with our objective selection criteria, to define the number of simulated 



- 15 - 



spectroscopic SN la confirmations. This number is divided by the number of input simulated 
SNe la to derive the yearly survey efficiency. 



4-2.1. Generating the Sample Population 
To simulate a realistic population of Type la SNe, we use the same LC templates and 



softw are used to determine photometric redshifts for the SNLS candidate SNe (ISullivan et al. 



20061 ). Figure [7] shows the canonical distributions of the parameters that characterize SN la 
LCs used in our efficiency simulations. The redshifts are chosen with a volume weighted 
uniform random number generator, to produce a redshift distribution over the range 0.2 < 
z < 0.6 that is uniform per unit volume as shown in Figure [TJi. We also calculated the run 
of dV(z), given the cosmological parameters from ^TJ and over-plotted this as a dashed line 
(with an offset of F = 0.005 for clarity) to show that our distribution is indeed constant per 
unit volume. The stretch values are selected using a Gaussian distribution centered on 1 with 
a width of a s = 0.1 (Figu re [7b). The intrinsic SN color is determined using the stretch-color 



relation from lKnop et al.l (120031 ). The host color excesses are chosen from the positive half of a 



Gaussian distribution, centered on 0.0 with a width of o~E(B-v) h — 0.2 ( Figure [7b). These ar e 



converted to host extinction assuming an extinction law with Ry = 3.1 (jCardelli et al.lll989l ). 



The peak magnitude offsets (after stretch correction) shown in Figu re [Tbl are chosen fro m a 



Gaussian distribution centered on with a width of <Jb max = 0.17 flHamuy et al.lll996l ). A 
uniform random number generator is used to pick the day of maximum for each simulated 
SNe la from a two- year-long interval that is centered on the middle of the survey range being 
simulated. This avoids problems with edge effects and produces an efficiency per year. We 
will address the systematic uncertainty due to differences between these distributions and 
the true distributions in §5.31 

In order to account for the possibility that a given SN can be missed because of tempo- 
rary localized losses of SN detectability in the Megacam array (see above, §4.1.21) . we assign 
a pseudo-pixel position to each simulated SN. This is done with a uniform random number 
generator that selects one of the 370 million Megacam pixels that are nominally available as 
the location of the SN. The number of real pixels available on a given epoch is calculated 
from the number of usable chips, derived during the fluxscale calculation. By choosing a 
random number out of 370 million, we are essentially assigning a probability that the SN 
will land on a region of the array that is temporarily unusable. If all chips are working, then 
the number of pixels available equals the nominal number and no SNe are lost. If a large 
number of chips are not working, the the number of pixels available is much less than the 
nominal number and a simulated SN has a higher probability of being missed. 
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4-2.2. Observing the Sample Population 

With the input sample population defined and LCs covering the simulated period gen- 
erated, we use the data describing the real SNLS survey epochs to observe the simulation. 
First, we use an a verage Milky Way extinction appropriate for the field being simulated 



(lAstier et al.ll2006l . Table 1). Then we use the epoch properties, equation [fl and Table E] 
to calculate visibility limits for each i' epoch. These visibility limits are used to define the 
signal-to-noise ratio, S/N, for each simulated SN observation using the following formula: 

S/N = 10.0 x ioI- .^-^ ^ 

with m e being the template magnitude of the simulated SN la in the epoch, and L e the epoch 
50% visibility limit. This formula assumes that an observation at the 50% visibility limit 
has a S/N of 10. The calculated S/N defines the width of a Gaussian noise distribution and 
a Gaussian random number generator is used to pick the noise offset for the observation. 

After the noise offsets are added to the observations, the resulting magnitudes in each 
epoch are then compared with the corresponding i' visibility limits, L e , and any magnitude 
that is brighter than its corresponding limit is considered an i' detection. We use a uniform 
random number generator to assign a real number ranging from 0.0 to 1.0 for each i'epoch. 
If this number exceeds 0.95, then the candidate is not detected in that epoch. This accounts 
for the 95% maximum recovery fraction observed in Figure [5K The shape of the recovery 
fraction at fainter magnitudes is already accounted for by using the 50% recovery magnitudes 
in the visibility limit calculation (see Figure \&&). To account for localized visibility losses, 
we calculate the number of pixels available on the epoch from the number of good CCDs 
available on that epoch. If the candidate was assigned a pseudo-pixel number larger than 
the number of good pixels on the epoch, then the candidate is not detected on that epoch. 

The restframe phases (relative to peak brightness) of all the relevant i' epochs are 
calculated for each simulated SN la using its given redshift. If a simulated SN la ends up 
with a detection in the restframe phase range from criterion 1 (see £J3]) then we evaluate it 
with respect to the remaining criteria. We calculate the restframe phase for each observation 
in the g', and r' epochs and then the remaining criteria are applied to decide if the simulated 
SN should be counted as a spectroscopically confirmed SN la. 

For the yearly efficiency, we keep track of the total number of SNe la that are simulated, 
since they are generated in yearly intervals. We also keep track of the number of SNe la that 
were simulated during the observing season for each field (from 158 to 211 days, see Tabled]). 
This allows us to compute our on-field detection efficiency and our on-field spectroscopic 
confirmation efficiency. 
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4-2.3. The Monte Carlo Survey Efficiency 

The resulting efficiencies for each field are presented in Tabled As we stated above, the 
statistical errors in these numbers are ~0.1%. We present the on-field i' detection efficiencies 
in column 2, which are all within 5% of 100%. This is expected considering the redshift range 
of our sample and the nominal %' frame limits. It also bolsters our spectroscopic completeness 
analysis by showing that our SN candidate list is not missing a significant population, in 
our redshift range. The on-field spectroscopic efficiency (column 3) averages close to 60%, 
reflecting the spectroscopic followup criteria applied to the detected SNe la. The yearly 
efficiency (column 4) averages close to 30% which reflects the half-year observing season for 
each field. 

We can compare Figure [3] with Table [7] as a consistency check. Starting with the on-field 
detection efficiencies (column 2), we notice that Dl has the lowest value. In Figure [3] we 
see that Dl has the largest variation in the visibility limits with some limits approaching 
i' = 20. Going to the spectroscopic on-field efficiency (column 3), we see that D2 has the 
lowest value. This is due to the large gap in the relatively short first season of D2. We also 
see that D4 has the highest on-field spectroscopic efficiency and the lowest scatter in the 
visibility limits in both seasons. In the last column of Table [7] we see that D3 has the highest 
yearly spectroscopic efficiency, due to the fact that D3 consistently has the longest seasons 
of any field (see also Tabled], column 7). D2 has the shortest season and consequently, has 
the lowest efficiency. 



5. Results 

We are now ready to apply our survey efficiencies to the culled, observed sample of 
SNe la and thereby derive a rate. The high detection efficiency of the survey from column 
2 of Table [7] illustrates that our sample for this study constitutes a volume limited sample 
as opposed to a magnitude limited sample. This means that we do not p roduce a predicte d 



redshift distribution to define our rate and average redshift as was done in lPain et al.1 (120021 ). 
for example. Instead, we apply our efficiency uniformly to our sample and our average 
redshift is the volume weighted average redshift in the range 0.2 < z < 0.6. We apply the 
appropriate efficiency to the sample of each field individually, propagating the Poisson errors 
of the field's sample, and then take an error-weighted average to derive our best estimate of 
the cosmic SN la rate averaged over our redshift range. We present the results from these 
calculations below. We also derive a rate per unit luminosity, present an analysis of our 
systematic errors, and compare our results with rates in the literature in this section. 
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5.1. SN Type la Rate Per Unit Comoving Volume 

We first need to calculate the true observed number of SNe la per year in each field. We 
then need to correct for the fact that at higher redshift, we are observing a shorter restframe 
interval due to time dilation. To derive the final volumetric rate we then calculate the total 
volume surveyed in each field and divide this out. We express these calculations with the 
following formula: 



rv = Nl ;' 2 [1 + (z)v\ {jj|^(0.6) - V(0. 

£yr ^SPEC 41/00 



•2)]}-\ (3) 



where Nj a /2 is the number of confirmed SNe la in the sample (see column 2 of Table [5]) 
divided by the number of seasons (2), e yr is the yearly spectroscopic efficiency from column 
4 of Table [3, Cspec is the spectroscopic completeness presented in column 6 of Table [5j 
1 + (z)v is the time dilation correction using the volume weighted average redshift over our 
redshift range, is the sky coverage in square degrees which is divided by 41253 (the total 
number of square degrees on the sky), and V(z) is the total volume of the universe out to 
the given redshift. These volumes are calculated using: 



V(z) = ~7T 



c f z dz' 



Ho Jo ^Vt m (1 + z'f + Q A 
with the parameters listed in £jl]and assuming a flat cosmology (flj. = 0). 



(4) 



The columns in Table [8] give the results at several stages in applying equation [3] along 
with some of the parameters used in the calculation. Column 2 presents the observed raw 
rate, r R Aw, calculated by simply dividing the average yearly sample for each field by the 
yearly spectroscopic efficiency, e yr . Column 3 shows r obs , the true observed yearly rate of 
SNe la in each field, which is the result of applying our spectroscopic completeness corrections 
{Cspec) to traw- Column 4 shows the results of accounting for time dilation by multiplying 
the observed rates by 1 + (z)y, where (z)y = 0.467 is the volume weighted average redshift. 
Column 5 lists the areal coverage for each field, after accounting for unusable regions of the 
survey images which include masked edge regions, and regions brighter than 20 mag arcsec -2 
in i' (see Figure [5H). The resulting survey volume between redshift 0.2 < z < 0.6 is reported 
in column 6, using equation H] which gives 1.035 x 10 6 Mpc 3 Deg -2 , using the cosmological 
parameters described in ^TJ Column 7 is the resulting volumetric rate. At each stage in the 
calculation the results are listed for each field and for a weighted average for the ensemble. 
Our derived rate per unit comoving volume, ry, is ry = 0.42 ± 0.06 x 10~ 4 yr _1 Mpc~ 3 
(statistical error only). 
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5.2. SN Type la Rate Per Unit Luminosity 



We use the galaxy L F derived from the first epoch data of the VIMOS-VLT Deep 
Survey (jllbert et al.ll2005l ) to calculate the B-band galaxy luminosity density. This recent 
survey derives the LF in the redshift range 0.2 < z < 0.6 from 2,178 galaxies selected at 
17.5 < Iab < 24.0. The Schechter parameters for the r est-frame B ban d LF are tabulated 
in their Table 1. We integrated the Schechter function (jSchechterl 119761 ) in the two redshift 
bins 0.2 < z < 0.4 and 0.4 < z < 0.6 and used a volume weighted average to get a luminosity 
density in the B-band of a B = 2.72 ± 0.48 x 1O 8 L , B Mpc" 3 . 

By using parameters derived from galaxies in our redshift range of interest, we do not 
need to evolve a local LF. As long as the slope of the faint end of the LF is well sampled 
and hence the a parameter is well determined, this produces an a ccurate luminosity density 
and hence an accurate SN la rate per unit luminosity. Figure 4 of lllbert et all (120051 ) shows 
that the LF in the highest redshift bin (0.4 < z < 0.6) is well sampled to ~3.5 magnitudes 
fainter than the 'knee' of the function. 



We now convert our volumetric rate into the commonly used luminosity specific unit 
called the "supernova unit" (SNu), the number of SNe per century per 10 10 solar luminosities 
in the rest-frame B band. Dividing the luminosity density in the rest-frame B band by our 
volumetric rate from Table M gives ri = 0.154+°;°^ SNu (statistical error only). 



5.3. Systematic Errors 

The rates for each field in Table [8] are all within one ustat of each other at each stage 
of the calculation of ry. This tells us that there are no statistically significant systematic 
errors associated with our individual treatment of the fields. In the subsequent analysis, we 
examine sources of systematic error that affect the survey in its entirety. We tabulate the 
values and sources of statistical and systematic errors in Table [9] for both ry and ri and 
describe each systematic error below. 



5.3.1. Spectroscopic Incompleteness 

We estimate the systematic error due to spectroscopic incompleteness by using our de- 
tailed examination of the SN candidates from §3.21 as tabulated in Table [51 column 4. Using 
the extremes of completeness for the ensemble (75% to 100%) as limits on this system- 
atic error, the spectroscopic incompleteness is responsible for a systematic error on ry of 
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(+0.03, -0.08) x 10~ 4 yr- 1 Mpc~ 3 and on r L of (+0.010, -0.031) SNu. 



5.3.2. Host Extinction 



For our canonical host extinction, we used a positive valued Gaussian E(B — V)h distri- 
bution wit h a width of o~r( R—y ) h = 0.2 (see Figure [7b) combined with an extinction law wit h 
Ry = 3.1 (jCardelli et al.lll989l ). We follow the procedure described in lSullivan et al.l (120061 ). 
with the exception that our host extinction is allowed to vary beyond E(B — V)h = 0.30. 
Systematics are introduced if our canonical distribution differs significantly from the real 
SN la host color excess distribution, or if there is evolution of dust properties such that 
the Ry = 3.1 model is significantly inaccurate. Preliminary results from submm surveys of 
SN la host galaxies out to redshift z = .5 show no significant evolution in the dust proper- 
ties when compared to hosts at z = 0.0 (IClements et al.l 120051 ) . We thus concentrate on the 
distribution of E(B — V)h as the major source of systematic errors in our redshift range. 

In an effort to quantify the systematic contribution of host extinction to an under- 
estimation of the SN la rate, we re-ran our Monte Carlo efficiency experiments setting 
E(B — V)h = 0.0 for each simulated SN. We analyzed the results of this experiment exactly 
as before (see §5.ip and derived a volumetric rate of ry = 0.38 + 0.05 x 10~ 4 yr -1 Mpc~ 3 . This 
rate is 10% (0.67<7stat) lower than the rate using our canonical distribution (see Table ED 
and quantifies the magnitude of the error possible, if host extinction is ignored. We will also 
use this zero dust rate to calculate rate corrections due to host extinction (see below). 

If we assume that our empirical host color excess distribution is biased by not including 
SNe in hosts with extreme E{B — V)h, and hence extreme Ay, then the systematic error 
on ry is positive only. In an attempt to quantify this err or, we compare our E (B — V)h 
distribution to models of SN la host extinction presented in iRiello fc Patatl (120051 . hereafter 
RP05). 



RP05 improve upon t he simple model of lHatano et al.l (119981 ) , motivated by the findings 
of Cappellaro et al. ( 1999 ) that the Hatano et al. ( 19981 ) model over-corrects the SN la rate 
in distant galaxies. RP05 use a more sophisticated model of dust distribution in SN la host 
galaxies and include the effects of varying the ratio of bulge-to-disk SNe la within the host. 
The resulting Ay distributions, binned by inclination, are strongly peaked at Ay = 0.0, have 
high extinction tails, and do not have a Gaussian shape (RP05, Figure 3). The smearing of 
the large fraction of objects with E(B — V)h ~ by photometric errors would produce a 
more Gaussian shape. 



Because the Ay distributions of RP05 and the Ab distributions of 



Hatano et al 



(11998 1 
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produce tails of objects with high extinction that extend beyond Gaussian wings, we per- 
formed two additional experiments using exponential distributions for E(B — V)h to simulate 
these tails. We generated these exponential distributions using a uniform random number 
generator to produce a set of random real numbers between and 1, which we will designate 
as 5i, and applied the following equation: 

E(B-V) h = -ln^/X E ( B -v) h , (5) 

where \e(b-v) h is the exponential distribution scale factor. The smaller the value of \E(B-v) h 
the larger the tail of the distribution. Figure [8] shows the two exponential distributions with 
^E{B-v) h — 5 and 3, along with the canonical distribution, converted to Ay using R v = 3.1 
and binned using the same bin size as RP05 (dAy = 0.1). If we compare these distributions 
with Figure 3 of RP05, we see that our canonical distribution is closest to the form of their 
model with an inclination range of 45° < i < 60°. Specifically, both our distributions have 
a maximum of Ay ~ 2.5. The exponential distributions are closer matches to their highest 
inclination bin, 75° < % < 90°, showing tails extending beyond Ay = 7.0 (although at very 
low probability). 

Using these distributions, we re-ran our Monte Carlo efficiency experiments and re- 
derived the volumetric rates to quantify the effect on our derived rate of missed SNe due to 
exponential tails in the host extinction distribution. For the \e(b-v) h — 5 distribution, we 
derived a rate of Ty = 0.44 ± 0.06 x 10~ 4 yr _1 Mpc -3 , which is only 5% higher than our 
canonical value. The \e{b-v) h = 3 case produced a rate of ry = 0.52 ± 0.07 x 10~ 4 yr _1 
Mpc" 3 , which is 24% or 1.67ctstat higher. This distribution is appropriate for spiral SN la 
hosts with high inclination, but will over-estimate the correction to rates from hosts with 
a range of inclinations and host morphologies. We, therefore, regard it as a measure of the 
upper limit on the statistical error due to host extinction. 

We can also compare the rate correction factors from RP05 with the correction factor 
resulting from the exponential dust distributions. Using our zero extinction experiment 
and the \e(b-v) h — 3 dust distribution, this factor is 1Z — 0.52/0.38 = 1.37. This value 
encompasses the factors reported in RP05 (see their §5) for their models with bulge-to-total 
SNe ratios of 0.0 < B/T < 0.5, which are given as 1.22 < TZb/t < 1.31. It also bounds 
their correction factors derived for dust models with Ry = 3.1, 4.0, and 5.0 which are given 
as 7Zr v = 1.27, 1.31, 1.34. Finally, this correction is not exceeded by the corrections derived 
from the RP05 models with total face-on optical depth, Ty = 0.5 and 1.0, which are given 
as TZ TV = 1.16,1.27. 

These comparisons demonstrate that it is reasonable to assume that we encompass the 
systematic errors due to host extinction if we use the \e(b-v) h — 3.0 exponential host 
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extinction distributions to define their upper limit. This results in a systematic error due to 
host extinction on ry of (+0.10) x 10~ 4 yr -1 Mpc~ 3 and on tl of (+0.037) SNu. 



5.3.3. Stretch 

We consider the effect of errors in the input stretch distribution on our derived rates. We 
re-ran our efficiency experiment doubling the width of the stretch distribution to a s = 0.2. 
This produced a rate of ry = 0.43 ± 0.06 x 10~ 4 yr _1 Mpc -3 , which is only 2% higher than 
the rate using a s = 0.1. The resulting systematic error due to stretch for ry is ±0.01 x 10~ 4 
yr -1 Mpc" 3 and for r L is ±0.040 SNu. 



5.3.4- Frame Limits 

Figure [6] shows the distribution of the 50% frame limits versus IQ e for a sample of %' 
images from the survey. The slope from Table El column 6 is steeper than one would expect 
from a simple analysis of the standard CCD S/N equation (e.g., Howell 1999) in the limit 
where the noise is dominated by the sky (i.e. at the frame limit). The expected slope is 
closer to 1.3, a factor of 1.7 lower than what we derived from our fake SN experiments. This 
could be the result of the coaddition process, perhaps because the co-alignment accuracy 
is sensitive to variations in IQ e . In order to account for a possible overestimation of the 
dependence of frame limit on IQ e (i.e. too large an a from equation [T]), we re-ran our 
efficiency experiment with a constant frame limit of i' = 24.5, which is shown in Figure [6] 
as the solid horizontal line. The resulting volumetric rate was ry = 0.39 ± 0.05 x 10~ 4 yr -1 
Mpc -3 . Using this value we estimate that an error in calculating the frame limits would 
introduce a systematic error on ry of (—0.03) x 10~ 4 yr -1 Mpc -3 and on of (—0.011) 
SNu. 



5.3.5. Host Offset 



One of the factors offered by iDahlen et al.l (120041 ) to account for the discrepancy be- 
tween ground-based rates near z = 0.5 and the delay-time models they present is the close 
proximity of SN candidates to host galaxy nuclei. If a candidate is too close to a bright 
host nucleus, they argue, it can be mis-classified as an AGN, or it might be passed up for 
spectroscopic followup because of the high level of host contamination. The results of our 
fake SN experiments (see §4. 1} Figure [6b) show that there is no such loss of SN sensitivity 
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close to the hosts of galaxies in the SNLS in the redshift range 0.2 < z < 0.6. Another way 
to look at these data is shown in Figure [9l which plots the percentage of fake SN missed 
as a function of host offset for two IQ e values from our fake SN experiments (see §4.ip . 
This further illustrates the lack of trend with host offset (see also Figure 2, Pain et al. 
2002). We specifically used a brightness weighted probability distribution for placing our 
fake SNe, shown in Figure [9] as the dot-dashed histogram, which preferentially places them 
in the brightest regions of a galaxy (i.e. near the center), so that we could detect any such 
problem. 



The study by iHowell et al.l (120001 ) also showed no significant loss of objects at small 
host offset when comparing a sample of 59 local SN la discovered with CCD detectors and 
a sample of 47 higher redshift (z > 0.3) CCD-discovered SNe. We conclude that this effect 
is not significant, at least out to z — 0.6, for the SNLS. 



5.3.6. Subluminous SN la Population 



Subluminous or SN 1991bg-like SNe la are another potential source of systematic error. 
Compared to the so-called normal SNe la these objects can have peak magnitudes up to 
two magnitudes fainter, exhibit different spectral features, have a different stretch color 
relation ship, yet they still obe y the stretch-luminosity relationship exhibited by the normal 
SNe la (jGarnavich et al.ll2004l ). As such, they would be useful on a Hubble diagram, however, 
none of our spectroscopically confirmed SNe la fall into the subluminous class. This has 
been confirmed in dependently by equivalent width measurements of our spectroscopically 
confirmed sample (lBronderll2006l ). We cannot assume, however, that the relative frequency 
of these objects decends to zero at higher redshifts. Their intrinsic faintness and differing 
color may produce a bias against these objects in our followup selection criteria. 

It is still likely that these subluminous SNe la are recent phenomena, since they are 
strongly associated with older stellar populations (lHoweilll200ll ). They have yet to be found 
in significant numbers at redshifts beyond z = 0.2, even though many CC SNe of similar 
peak magnitude have been found. It is thus reasonable to assume that the current fraction 
of subluminous SNe la forms an upper limit on the f raction at higher redshifts. The best 
estimate of the current fraction is from lLi et al.l (120011 ). who derive a fraction of 16 ± 6% for 
the subluminious class using a volume-limited sample. Rather than attempt to evolve this 
number to the redshift range of interest for this study, we take the conservative approach 
and use it to calculate an upper limit on our systematic error due to missing the subluminous 
SNe la, assuming their relative fraction has no evolution. This calculation yields a systematic 
error on r v of (+0.08) x 10~ 4 yr" 1 Mpc" 3 and on r L of (+0.029) SNu. 
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5.4. Comparison with Rates in the Literature 



Figure [TO] shows the same rates and SFH as Figure [fl but now with our rate included 
as the filled square. At z ~ 0.45 all four observed rates are in statistical agreement. 



The result from lPain et al.l (120021 ) at z — 0.55 of ry 



0.525lg5S(stat) 



^:^(syst)xio- 4 



yr Mpc is h igher than our rate, but still in statistical agreement. We must consider 
that iPain et al.l (120021 ) do not account for host extinction in their rate derivation (see their 
§6.8). There are two possible explanations. Either host extinction has a small effect in 
calculating SN la rates in this redshift range, or the lack of host extinct ion correct i on wa s 
compensated for by an equal amount of contamination in the result from lPain et al.l (120021 ) . 
If we calculate a correction for host extinction from our Monte Carlo experiment where we 
set all SNe to have E(B — V)h = 0.0, we can estimate how their rate would change if they 
had accounted for host extinction using our method. This correction is 1Z = 0.42/0.38 = 1.11 
when compare d to our canonica l host extinction results. Applying this 11% correction to 
the value from lPain et al.l (120021 ) produces a rate of ry = 0.58 x 10~ 4 yr -1 Mpc -3 , which is 
only 0.5<t higher than their original value. We, therefore, conclude that the rates derived in 
this range are not significantly affected by host extinction. 



5.4-1- Contamination 



Even though all t he rat es near z = 0.5 are in statistical agreement, the rate at z — 0.55 
from iBarris fc Tonryl (120061 ) is within our redshift range and is nearly five times our value 
(> 4cr greater). The largest disagreement between published rates in Figure ITUl ( and i n the 
literature as far as we know) is that betw een the rates at z = 0. 5 5 of IPain et al.l (120021 ) and 



B arris fc To nrv (2006). We point out that lBarris fc Tonryl (120061) is a re-analysis of the data 



from [Tonrv et al.l (120031). wh ich reported a rate that agrees with IPain et al.l (120021) . The rat e 
from IBarris fc Tonryl (120061 ) is a factor of 3.9 higher than the rate from IPain et al.l (120021 ). 
If we add the errors on these rates in quadrature, this amounts to a 3.8<r difference. 

We have shown that host extinction cannot be the sole e xplanation f or the se discrepan- 
cies. Our estimate for the host extinction correction factor for lPain et al.l (120021 ) is 72. = 1.11, 
which would not be enough to resolve it. A host ext i nction correction factor of 1Z = 2.61 
would be required to bring the rate from IPain et al.l (120021 ) just in to statistical agreemen t 



with the rate of ry = 2.04±0.38 x 10~ 4 yr -1 Mpc" 3 at z = 0.55 from IBarris fc Tonrvl (120061 ). 



This correction is larger than the model from RP05 for hosts with total face-on optical depth 
of Ty = 10 which gives TZ(ry) = 2.35, the largest correction listed in RP05. An even larger 
correction would be required to bring our rate into statistical agreement with the rate from 
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B arris &: Tonryl (120061 ) at z — 0.55. We must await deeper submm studies to see if correc- 
tious that large are r easonable for SN la hosts out to redshift z = 0.6. Indications from 
Clements et al.l (120051 ) are that this does not describe SN la hosts out to redshift z = 0.5. 



We maintai n that contamin a tion i s the largest source of systematic error in SN la rates 
beyond z = 0.5. iBarris fc Tonryl (120061 ) us ed LCs generated f rom relatively sparsely sampled 
(At = 2-3 weeks) RIZ filter photometry (IBarris et al.l 120041 ) co mbined with a traini ng set 
of 23 spectrally identified SN la to verify their SN la sample. iDahlen et al.l (120041 ) used 
low-resolution grism spectroscopy i n combination with p hotometric methods to identify the 
majority of their candidate SNe la. IStrolger et al.l (120041 ) sta te that luminous SNe Ib/ c can 



occupy nearly the same magnitude- color space as SNe la. Uohnson fc Crottsl (120061) also 
conclu de that SNe Ib/c are the biggest challenge in phototyping SNe la. Istrolger et al. 
(J200J) point out that the bright SNe Ib/c make up only ~ 20% of all S Ne Ib/c and that 



SNe I b/c make up only ~ 30% of all core-collapse (CC) SNe according to ICappellaro et al. 



( 119991 ) . One obvious caution is the fact that these ratios are based on a small sample (< 15) 
from the local universe. Star-formation increases with redshift and it is plausible that the 
relative frequency of SNe Ib/c may increase as well. 

We also assert that CC SNe at lower redshifts can masquerad e as SNe la at highe r 



redshifts. We support this assertion by pointing out Figure 9 in ISullivan et al.l (120061 ). 



This figure plots photometrically determined redshifts (using a SN la template) against 
spectroscopically determined redshifts for SNe la and several types of CC SNe. The redshifts 
of the CC SNe are systematically over-estimated by as much as Az = +0.5. Contamination 
of this kind is not addressed by using the typical m agnitude difference between CC SNe and 
SNe la to cull out CC SNe (IRichardson et al.l 120021 ). since the CC SNe appear at the wrong 
redshift. This problem is strongly mitigated when host redshifts are available to cross-check 
SN photo-redshifts, as long as the host redshifts are reasonably accurate. 

Another problem for photometric typing is reddening. As mentioned above, SNe la are 
distinguished from CC SNe in the early part of their LCs by having redder colors. Thus, 
a highly reddened CC SNe can appear to be a SN la and weeding these objects out of a 
SN la sample requires good epoch coverage of the later epochs. The SNLS has the benefit 
of 4 filter photometry which helps distinguish even highly reddened CC SNe early on. Even 
with this advantage, ~ 10% of our c andidates promoted for spectroscopic followup turn out 
to be an identifiable type of CC SN (IHowell et al.l 120051 ) . 



The diversity of CC SNe, as compared to SNe l a, is another challenge for photomet- 
ric identification of SNe. Neither IDahlen et al.l (120041 ) nor IBarris fc Tonryl (120061 ) include a 
large database of spectrally identified CC SNe light curves in their training or test data sets. 
Until photometric methods can prove themselves convincingly against the full diversity of 
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CC SNe, spectroscopy is the most reliable way to identify SNe. The pay-off for develop- 
ing an accurate classifier based only on photomet ry is huge, however , given the expense of 
obtaining spectroscopy of SNe at high redshifts. IJohnson fa Crottsl (120061 ) point out that 
having good spectral-energy distribution coverage and having dense time-sampling will im- 
prove the accuracy of this method, a state ment that agrees with our experience with the 
SNLS JSullivan et al.lbood : buv et al.ll2005h . 



We have emphasized the importance of verifying the majority of the SN la sample with 
spectroscopy in all of our figures comparing rates and SFH by plotting all the rates from 
samples that fulfill this criterion as filled symbols. Assuming these surveys have carefully 
characterized their spectroscopic completeness, the trend they display is the one that should 
be compared to SN la production models. 



If we look in detail at the lDahlen et al.l (120041 ) sample, we point out that the two redshift 
bins in which the m ajority of objects are spectroscopically confirmed (the 'Gold' objects, see 
Strolger et al.l 120041 ) at z = 0.5 (2 out of 3) and z = 1.2 (5 out of 6) are included in our 
comparison of SFH and SN la rate evolution (see Figure [TU]) . The other two bins have only 
50% of their sample in the 'Gold' category: 7 out of 14 at z = 0.8, and 1 out of 2 at z = 1.6, 
and are therefore not included in this comparison. 



6. Discussion 

Our rest-frame SN la rate per unit comoving volume in the redshift range 0.2 < z < 0.6 
((z)v — 0.47) and using the cosmological parameters Q m = 0.3, = 0.7 and H = 70 km 
s _1 is 

r v ((z) v = 0.47) = 0A2t° o s o l{syst) ± 0.06(stat) x l^yr^Mpc^ , (6) 

and we also report our SN la rate in SNu, for comparison with previously determined rates 
at lower redshift: 

r L ({z)v = 0.47) = OAMt Z(syst)t° Z(stat)SNu. (7) 



6.1. Comparison with Star Formation History 

The place to begin investigating the relationship between SFH and SN la production is 
where the systematic errors are minimized. Volumetric rates are the most appropriate to use 
in exploring this relationship. The SNu, defined using a blue host luminosity, is less ideal 
especially for SNe la for a number of reasons. Galaxy luminosity evolution makes interpreting 
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trends in SNu with redshift difficult. Also, using a blue luminosity is not good for SNe la, 
since they also appear in galaxies with older and redder populations than CC SNe. 

The region of minimal systematic uncertainty in the volumetric rate evolution is de- 
lineated by the trade-off between survey sensitivity and volume sampling. At low redshifts 
most searches are galaxy-targeted, requiring the conversion of the luminosity specific rate 
(SNu) to a volumetric rate through the local galaxy luminosity function. The volume sam- 
pled is low, increasing the influence of cosmic variance on the derived volumetric rate and 
hence, increasing the systematic errors. At high redshifts, survey sensitivity dominates the 
systematics since high redshift SNe are close to the detection limits, must be spectrally con- 
firmed with lower S/N spectra, or photometrically identified with lower S/N LCs, and have 
projected distances that are fewer pixels from host nuclei. 

For the SNLS, the redshift range 0.2 < z < 0.6 is low enough to reduce systematic 
errors associated with detection limits, spectral confirmation, and host offset. It also sam- 
ples 4 large, independent volumes of the universe minimizing the effects of cosmic variance 
on the errors. In our subsequent analysis, we examine two currently popular models of 
SN la production and compare them to our SN la rate at (z)v = 0.467, and to the other 
spectroscopically confirmed rates from the literature. 



6.1.1. The Two-Component Model 



This recent model, first put forth by iMannucci et al.l (20051 ) and applied to a sample of 
rates from the literature by Scannapieco fe Bildsten " ([20051 ) . proposes a delay function with 



two components. One component, called the prompt component, tracks SFH with a fairly 
short delay time (< 1 Gyr). The other component, called the extended component, is pro- 
portional to total stellar mass and has a much longer delay time. This model arose as a way to 
account for the high SN la rate in actively star-forming galaxies, relative to less active galax- 



(Oemler & Tinslev 


19791; 


van den Bereh 


l99o|; 


Cappellaro et al. 


1999: } 


vlannucci et al. 


2005; 


Sullivan et al. 


2006 


)• 


Gal- Yam & Maoz 


(2004 


) and 


Maoz & Gal- Yam 


( 2004J) observed the 



SN la rates in galaxy clusters and indicated the possibility that a lo ng delay-time may be in 



consis tent with their observations, given the cluster Fe abundances. IScannapieco fe Bildsten 



( 120051 ) demonstrate that the Fe content of the gas in galaxy clusters can be explained by 
the prompt component of the two-component model. They also demonstrate that the two- 
component model reproduces t he observed stellar [O /Fe] abundance ratios within the Galaxy 



( IScannapieco fe Bildsten! 120051 . Figure 3). 
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The two-component model is described by equation 1 from IScannapieco fc Bildsten 



(120051 ). which gives the relationship between the cosmic star formation rate, the cosmic 
stellar mass and the volumetric SN la rate as a function of time. We re-write their equation 
in a more general form here: 



r v (t) = AM*(t) + BM*(t), 



(8) 



which gives A in terms of SNe la per year per unit mass and B in terms of SNe la per year 
per unit star formation. The component scaled by A is the extended component, while B 
scales the prompt component. 

By comparing this model to our rate, and other spectroscopically confirmed rates from 
the literature, we can estimate the relative contributions of the extended and the prompt 
components. Figure ITUl already places an upper limit of B < 1CT 3 yr~ 1 (M Q yr^ 1 )^ 1 because 
the SFH curve normalized to go through our rate represents a pure prompt component model 
(A = 0.0). 

Figure [11] shows a non-linear least squares fit of the spectroscopically confirmed SN la 
rates from this study (filled square) and the lit erature (filled circles) to the rates predicted by 
the two-component model using the SFH from lHopkins fc Beacoml (120061 ). This fit produces 
a reduced \ 2 goodness-of-fit statistic of xl — 0.361 using the published error bars, produced 
by adding the statistical and systematic errors in quadrature. With 7 degrees of freedom (9 
data points minus two parameters), this corresponds to a probability of 93% that the null 
hypothesis is correct, i.e., that the data represent a random sampling from the parent distri- 
bution described by the fit. The resulting fit parameters are A 
and B = 8.0 ± 2.6 x 10~ 4 yr~ 1 {M Q yr' 1 )- 1 (see Figure HI]). 



1.4 ± 1.0 x 10~ 14 yr^M' 1 



This is t he first time this model has b een fit directly to volumetric rate data from 
the literature. IScannapieco fc Bildsten! (120051 ) normalized each component separately (using 
a different SFH) and then compared the resulting rate evolution to the observed rates. 
For the extended (A) component, they used the rate per unit mass for E/S0 galaxies from 



Mannucci et al.l (120051 . Table 2) which gives a value ofA = 3.8+^x 10" 14 yr^M^ 1 , corrected 



to o ur cosmology. An alte rnative value for the A component can be derived from Table 
3 in iMannucci et al.l (120051 ) which gives the rate in bins of B — K color, independent of 
morphology. The reddest bin, having B — K > 4.1 gives A = 2.4^^ x 10~ 14 yr -1 !^^ 1 
(again adjusted for our cosmology), which is more consistent with the value from our fit. A 
possibly more important difference stems from the definition of mass. Our mass is derived 
from the integration of the SFH from high redshift to the epoch in question and , there fore, 
includes the mass from stars that have died. The mass used in IMannucci et al.l (120051 ) was 
derived from the i^-band luminosity of individual galaxies and is more representative of the 
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mass currently in stars. Since our method tends to over-estimate the mass, our A value is 
correspondingly lower. 



For the prompt (B) component, IScannapieco fc Bildstenl (120051 ) first normalize the CC 
rate to the SFH and then use an assumed CC/SN la ratio, a method that they admit is highly 
uncertain. This produces a prompt component with B = 23 ± 10 x 10 -4 yr~ 1 (M & yr -1 ) -1 , 
which is 1.45a higher than our value. They also mention an alternate normalization using 
the SN la rate in actively star-forming galaxies as indicated by having B — K < 2.6. This 



method produces the value of B 
our value. 



IOI5 x 10 4 yr 1 (M yr x ) 1 , which is consistent with 



We are comparing a met hod that normalizes the A and B components separately 
(IScannapieco fc Bildstenl 120051 ) . with a method which directly ties the component values 
to the rate evolution, as delineated by the rates derived from spectrally confirmed samples. 
Given the differences in SFH, mass definition, and method, the level of agreement is encour- 
aging for this model. Caveats remain, however, including systematics in the lower redshift 
rates due to cosmic variance and systematics in the SFH. Further tests of this model will 
come as the SN la rate evolution is more precisely measured (at low and high redshifts) and 
as the SN la rate per unit mass an d per unit star-form ation is measured more accurately for 
a larger set of individual galaxies (jSullivan et al.ll2006l ). 



The success of this model implies that it is reasonable to des cribe SN la productio n in 
terms of two populations with two different delay times (see also iMannucci et al.ll2006l ). If 
these two populations represent two separate channels for SN la production, they may also 



exhib it different intrinsic properties. The SNLS is carefully examining this ([Sullivan et al. 



20061 ) to avoid biases in our cosmological parameters. At higher redshifts the component tied 
to SFH will tend to dominate, while lower redshift samples will contain more of the extended 
component SNe la. Cosmological parameters determined with SNe la spanning a large range 
of redshifts may be subject to systematics, if unaccounted-for intrinsic differences exist. 
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6.1.2. Gaussian Delay Time Model 
Figure [121 shows a comparison of observed rates with the delay time model delineated 



Strolger et al.l (120041 ). which is plotted as a dashed line. This model convolves the SFH 
with a Gaussian delay-time distribution with a characteristic delay time, r, and a width 
that is some frac tion of the delay time: a = 0.2r in this case. We updated the SFH model, 
using the fit f r omlHopkins fc Beacoml (120061 ). We find that the delay time model still fits the 
Dahlen et al. torn data with a delay time of r = 3 Gyr, which is statistically consistent 
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with the most like ly value from IStrolger et al.l (12004 ). In contrast, iMannucci et al.l (120061 ) 
used the data from lDahlen et al.l (120041 ) and combined with host galaxy colors and host radio 
properties found a bi-mod al delay distribution to be more consistent. We find that the single 
Gaussian model fit to the iDahlen et al.l (120041 ) d ata consistently o ver-p redicts the rates near 
and below z = 0.5. In particular, the rate from iBarris fc Tonryl (120061 ) at z = 0.25 is more 
than la below this model. 



We also show a Gaussian delay time model normalized to our rate in Figure d2] as the 
solid line. It is statistically consistent with the spectrally confirmed SN la rates, except the 
highest one at z — 1.2. It has the following parameters: r = 4.0 Gyr and a = 0.7r. This 
model predicts a very low SN la rate at higher redshifts in contrast to the photometrically 
typed rates near a redshift of z ~ 0.7. 

While the observations can be fit with this model, the favoured delay times tend to 
be longer than 3 Gyr. This is inconsistent with the finding that th e SN la rate is much 
higher in galaxies with recent star formation (Qemle r fc Tinslevlll979l ; Ivan den Berghlll990 



Cappellaro et al. 1999 ; Mannucci et al. 2005) ; 



ferent set of observed rates, IMannucci et al.l (120061 ) find that a single Gaussian delay time 



Sullivan et al. 20061 ) . Although we use a dif- 



does not fit observed rates as a function of redshift, host color and radio loudness as well as 
a bi-modal delay-time distribution. If there is a strong correlation in the SN la rate with 
host galaxy star formation rate, then the overall SN la rate evolution should track the SFH 
reasonably closely, especially at higher redshifts. 

The real test of the Gaussian delay model will come with rates beyond z = 1.4 where the 
predicted down-turn compared to SFH w ill become pronounced. Our current best estimate 
in this range is from IDahlen et al.l (120041 ) and is based on a sam ple of two SNe, only one of 
which was a member of the 'Gold' set from lStrolger et al.l (120041 ) . A larger, spectroscopically 
confirmed sample fr om a much deeper survey is needed to confirm or refute this result from 
Dahlen et al.l (120041 ) and therby either support or discredit the Gaussian delay-time model. 



This model faces another challenge. If there is a large component of the SN la rate that 
is closely tied to the SFH, then at higher redshifts the majority of SNe la will arise closer to 
star-forming regions in their hosts. This implies that the effect of host extinction (dust) on 
the sensitivity of SN la surveys will grow with redsh ift. It must be shown that the downturn 
in the SN la rate measured by IDahlen et al.l (120041 ) and predicted by this model is not the 
result of these factors. To do this requires the use of updated host extinction models such 
as RP05 combined with dust evolution models, derived from the deepest IR and submm 
surveys. 
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7. Summary 

We have produced the most accurate SN la rate to date by using a spectroscopically 
confirmed sample and detection efficiencies derived from a well characterized survey. We 
investigated known sources of systematic errors using recent models of host extinction from 
RP05, and fake SN experiments to test host contamination losses. Our derived volumetric 
rate from a culled sample of 58 SNe la in the redshift range 0.2 < z < 0.6 is ry = 0.42l ; g 
(syst) ±0.06 (stat) xlO -4 yr _1 Mpc -3 . We conclude from our experiments and a comparison 
of other rates in the literature, that contamination may be the largest source of systematic 
error for rates up to redshift z — 1, in particular, for those rates based on samples that are 
photometrically typed. 



Using the recent SFH fit fr om iHopkins fc Beacoml (120061) . we compare our rate with 



the two-component model from IScannapieco fc Bildstenl (120051 ) and place an upper limit 



on the contribution from the component of SN la production that is closely tied to star 
formation of B < 10~ 3 yr~ 1 (M & yr^ 1 )^ 1 . By fitting this model to our rate and the spectrally 
confirmed rates in the literature we make an estimate of both components directly and find 
A = 1A± 0.9 x 10~ 14 yr- l M~ l and B = 8.0 ± 2.4 x 10" 4 yr~ 1 {M Q yr' 1 )- 1 , with the caveat 
that our mass definition is an over-estimate (is the integral of SFH and, therefore, includes 
dead stars). 
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Fig. 1. — Observed SN la rate versus redshift. The filled circles are SN la rates derived 
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Fig. 2. — Nighly averaged photometry for the spectroscopically confirmed sample of SNLS 
SNe la. Each filled circle represents the nightly average magnitude for the g' (green), r' 
(orange), i' (red), and z' (purple) filters on a normalized AB magnitude scale. The template 
fits for each filter are indicated by the solid lines of the corresponding color. The designation 
from Table [2] (minus the SNLS prefix) is indicated in the upper left corner and the spec- 
troscopic redshift is indicated in the upper right corner of each panel. If the SN was culled 
from the sample, this is indicated by the word 'Rej' under the designation. The remaining 
SNe are presented in the following pages. 
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Fig. 3. — Julian Date versus magnitude limit for the i' epochs in each of the four deep fields 
used for SN la detection in this study. The calendar year transitions are indicated by the 
vertical dashed lines. The i' magnitude limits plotted are described in §4.11 
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Fig. 4. — Nighly averaged photometry (as in Figure W) for the unconfirmed SNLS SNe la 
from Table HI The designation is in the upper left corner and the photometric redshift is in 
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designation for each SN. The first six objects are presented here. The full set is available 
from the online version of the article. 



-54- 



Restframe Day 
-20 -10 10 20 



20 



CD 

3 21 



CC 



22 



CD 
N 



« 23 



24 



-40 -20 



30 



05D3lx 
X 2 : 3.lf 


^z ph :0.58 


i / ~» 






\ 


N / W 





20 40 



Restframe Day 
-20 -10 10 20 



CD 
"O 

"c 
O) 
ctf 

E 

"O 
CD 
N 

Ta 
E 



20 
20 
21 
21 
22 
22 



03D4bx 
% 2 : 2Arf*-% 



z ph : 0.54; 




-20 20 40 



-20 -10 10 20 -20 -10 10 20 30 




-40 -20 20 40 -20 20 40 




4- 



continued. 



-55 - 



Restframe Day 
-20 -10 10 20 30 



04D1fh_^ 
X 2 : 5.a?% 


z ph :0.56: 


// /-"' ^ 


fjj 




'■ I 1 /- r ' 

1 i / -9 


\ < 





CD 
"O 

'c 

O) 
ctf 

E 

"O 
CD 
N 



CO 

E 
i— 
o 



Restframe Day 
-20 -10 10 20 30 



20 



21 - 



£ 22 



23 



04D1qn 
% 2 : 6.4 


^ z ph : 0.52; 


k if 


%\ • 


* // 

\h 1L 

M 5 IP 


\ V i 

A i : 



-20 



20 



40 



-40 -20 



20 40 




-20 -10 10 20 30 



-20 -10 



10 20 30 



L 05D3bz»- 

;% 2 :5^* 


Z ph : 


0.52^ 






i 


i I / 
J 


\ ^ 


V \ 









-20 20 40 
Observed Day 



20 

CD 

1 20 

"c 
CO 

i 21 
■a 

CD 

I 21 

E 
i— 
o 

z 22- 



05D3ch 
X 2 : 5.5 



z ph : 0.35 




-20 -10 10 20 30 
Observed Day 



4- 



continued. 




4. — continued. 



-57- 



Q 
LU 
I— 

o 

LU 
I— 
LU 
Q 



100 
80 
60 
40 
20 














i 










1 
r 

/ 

# 


0.69" 




/ 

/ 


1.06" : 




/ 





24 23 22 21 
SN MAG (i') 



100 
80 
60 
40 
20 






i i 
fa d i 




b): 


: i' < 50% LIMIT 









2 4 6 8 10 
HOST OFFSET (ASEC) 



Q 

LU 
I— 

o 

LU 
I— 
LU 
Q 



100 






100 






80 


_ i 


" 1 P | * _ | ' 1 — 1 


80 


|- . 




60 




c). 


60 




d): 


40 


. "L 




40 






20 


i' < 50% LIMIT 


20 


i' < 50% LIMIT 
















■ 




-2 


-10 12 




28 26 24 22 20 


18 






HOST - SN MAG 




BKG (MAG/ASEC A 2) 



Fig. 5. — Recovery percentage after human review for the two i' images having IQ e = 0'.'69 
(solid line) and IQ e = 1'.'06 (dashed line): a) shows the total recovered percentage as a 
function of fake SN i' magnitude with the 50% recovery limits for each IQ shown as the 
vertical lines, b) - d) show the recovery fraction versus various parameters for the fake SN 
above the 50% recovery fractions: b ) shows the recovered percentage as a function of host 
offset in arcseconds c) shows the total recovered percentage as a function of host minus 
fake SN %' magnitude, and d) shows the recovered percentage as a function of %' background 
measured in i' magnitudes per square arcsecond. 
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Fig. 6. — 50% detection limit as a function of IQ e for a range of SNLS i' images with a range 
of IQ e . The IQ e distribution for all i' images in the survey sample time range is shown as 
the dashed histogram. Automated detection limits for field Dl are indicated by the open 
diamonds, and for field D4 by the open squares. The 50% recovery limits after human review 
for two Dl images at IQ e = 0'.'69 and 1'.'06 are indicated by the asterisks. A linear fit to the 
human review limits is indicated by the bold solid line. A constant frame limit, independent 
of IQ e , is indicated by the horizontal thin line at i' = 24.5. The human review fit and the 
constant limit at i' = 24.5 encompass all the limits shown and define the range used in our 
calculation of the systematic errors associated with SN la detection in our simulations (see 
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Fig. 7. — Canonical distributions of properties of the simulated SNe la used in the Monte 
Carlo efficiency experiments: a) shows the volume- weighted redshift distribution as the 
solid histogram and the run of dV(z) as the dotted line, b) shows the Gaussian stretch 
distribution with a = 0.1, c) shows the positive- valued Gaussian host extinction distribution 
with aE(B-v) h = 0.2, and d) shows the Gaussian ABmax distribution with &b max = 0-17 
with a Gaussian fit overplotted as a thin solid line. The fitted a of the Gaus s ian is annotated 
on the plot and matches the distribution of real SNe la from lHamuy et al.l (119961 ). 
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Fig. 8. — Distributions of total V-band extinction, Ay, for three models of SN la host 
extinction. The solid line represents the positive- valued Gaussian with a width of <Je(b-v) h — 
0.2 that was used for the canonical Monte Carlo efficiency experiments. The dashed line 
represents an exponential distribution of E(B— V)h with a scale parameter of \e(b-v) h — 5-0. 
The dotted line shows an exponential distribution with a scal e parameter of \F,m-y ) h = 3.0. 
When comparing these distributions to those in Figure 3 of iRiello fc Patatl (120051 ) . we see 
that our canonical host extinction model is appropriate for an intermediate host inclination 
model (45° < i < 60°), while the exponential distributions are closer to the extreme host 
inclination model (75° < i < 90°). 
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Fig. 9. — Percent missed versus host offset in arcseconds for IQ e = 0"69 (solid line) and 
IQ e = 1'.'06 (dashed line) from the fake SN experiments described in §4.11 The error bars 
are the Poisson errors in each bin. The histogram of the input host offsets (divided by 10) is 
plotted as the dot-dashed line. All data were binned with 0"25 bins. The percentage missed 
does not turn upwards at low host offset as would be expected if there was a loss of SN 
visibility near the host galaxy nuclei. 
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Fig. 10. — Same as Figure [fl but with the ra te from this study p l otted as a filled square. 
Because the renormalization of the SFH from Ifiopkins fc Beacoml (120061 ) using a factor of 
10 3 fits our rate, we can immediately place an upper limit on any component of SN la 
rate production that is tied directly to star formation of t S fh ^ 1 SN la / 10 3 M Q (B < 
10 -3 yr ~ l (M Q yr' 1 )- 1 , see 36~TTj) . 
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Fig. 11. — Fit of observed, spectrally confirmed SN la rates (see Figures IT1 and [TUl) to the two 
component model from equation [8] shown as the solid line, with the one-sigma errors defining 
the shaded region. The extended (A) component is proport ional to the current mass density, 
defined by integrating the SFH of lHopkins fc Beacoml (120061 ). and is shown by the dotted line, 
while the the prompt (B) component is proportional to the instantaneous SFH and is shown 
by the dashed line. The non-linear least-squares fit to the spectroscopically confirmed rates 
has xl = 0.361, and produces an extended component with A = 1.4 ± 1.0 x 10~ 14 yr^M^ 1 
and a prompt component with B = 8.0 ± 2.6 x 10~ 4 yr~ 1 (M Q yr^ 1 )^ 1 . The errors quoted 
are statistical only and do not include systematics due to errors in the SFH or the mass 
definition. 
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Fig. 12. — C o mpar ison of observed SN la ra tes with the delay-t i me m odel reported in 



Strolger et al.l (120041 ) . but using the SFH from iHopkins fc Beacoml (120061 ). The filled and 



open symbols are t he same as in Figure dU The dashed line is a fit to the data from 
Dahlen et al.l (120041 ) with parameters of r = 3 Gyr, and a = 0.2r. The solid line is a fit that 
is no rmalized at our rate having parameters of r = 4 Gyr, and a = 0.7r. The fit to the data 
from iDahlen et al.l (120041 ) systematically over-predicts the rates at lower redshift. Both fits 
predict low SN la rates beyond z = 1.5. 
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Table 1. SNLS Observations 



Field 


Position(J2000) 


Season 




Start Date 




End Date 


Span 


RA 


Dec 


(MJD) 


(YYYY-MM-DD) 


(MJD) 


(YYYY-MM-DD) 


(days) 


Dl 


02:26:00.00 


-04:30:00.0 


1 


52852 


2003-08-01 


53026 


2004-01-22 


174 








2 


53207 


2004-07-21 


53390 


2005-01-20 


183 


D2 


10:00:28.60 


+02:12:21.0 


1 


52993 


2003-12-20 


53151 


2004-05-26 


158 








2 


53328 


2004-11-19 


53503 


2005-05-13 


175 


D3 


14:19:28.01 


+52:40:41.0 


1 


53017 


2004-01-13 


53228 


2004-08-11 


211 








2 


53386 


2005-01-16 


53586 


2005-08-04 


200 


D4 


22:15:31.67 


-17:44:05.0 


1 


52795 


2003-06-05 


52964 


2003-11-21 


169 








2 


53173 


2004-06-17 


53350 


2004-12-11 


177 
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Table 2. Spectroscopically Confirmed Type la SNe: 0.2 < z < 0.6 



Name QJ2000.0 5j2000.0 zspec Discovery (MJD) Culling Status a Refs b 



SNLS-03Dlar 


02:27:14.680 


-04:19:05.05 


0.41 


52901 






SNLS-03Dlau 


02:24:10.380 


-04:02:14.96 


0.50 


52901 




A06 


SNLS-03Dlaw 


02:24:14.780 


-04:31:01.61 


0.58 


52901 




A06 


SNLS-03Dlax 


02:24:23.320 


-04:43:14.41 


0.50 


52901 




A06,H05 


SNLS-03Dlbp 


02:26:37.720 


-04:50:19.34 


0.35 


52913 




A06,B06 


SNLS-03Dldj 


02:26:19.082 


-04:07:09.38 


0.40 


52962 


1-3,5 


A06 


SNLS-03Dlfb 


02:27:12.855 


-04:07:16.40 


0.50 


52991 




A06,B06 


SNLS-03Dlfc 


02:25:43.602 


-04:08:38.77 


0.33 


52991 




A06,B06 


SNLS-03Dlgt 


02:24:56.012 


-04:07:37.08 


0.55 


53000 


2,4 


A06,B06 


SNLS-04Dlag 


02:24:41.108 


-04:17:19.69 


0.56 


53019 


5 


A06,B06 


SNLS-04Dlak 


02:27:33.399 


-04:19:38.73 


0.53 


53019 




A06.B06 


SNLS-04Dldc 


02:26:18.477 


-04:18:43.28 


0.21 


53228 


2 




SNLS-04Dlhd 


02:26:08.850 


-04:06:35.22 


0.37 


53254 




H05 


SNLS-04Dlhx 


02:24:42.485 


-04:47:25.38 


0.56 


53254 






SNLS-04Dljg 


02:26:12.567 


-04:08:05.34 


0.58 


53267 






SNLS-04Dlkj 


02:27:52.669 


-04:10:49.29 


0.58 


53283 






SNLS-04Dloh 


02:25:02.372 


-04:14:10.52 


0.59 


53294 






SNLS-04Dlpg 


02:27:04.162 


-04:10:31.35 


0.51 


53312 






SNLS-04Dlrh 


02:27:47.160 


-04:15:13.60 


0.43 


53344 






SNLS-04Dlsa 


02:27:56.161 


-04:10:34.31 


0.59 


53351 






SNLS-04D2ac 


10:00:18.924 


+02:41:21.45 


0.35 


53022 


5 


A06,B06 


SNLS-04D2bt 


09:59:32.739 


+02:14:53.22 


0.22 


53081 


1-4 


A06,B06 


SNLS-04D2cf 


10:01:56.048 


+01:52:45.90 


0.37 


53081 


1-4 


A06,B06 


SNLS-04D2cw 


10:01:22.821 


+02:11:55.66 


0.57 


53081 


2-5 


A06,B06 


SNLS-04D2fp 


09:59:28.183 


+02:19:15.20 


0.41 


53094 




A06.B06 


SNLS-04D2fs 


10:00:22.110 


+01:45:55.64 


0.36 


53094 




A06,B06 


SNLS-04D2gb 


10:02:22.712 


+01:53:39.16 


0.43 


53094 




A06 


SNLS-04D2gc 


10:01:39.267 


+01:52:59.52 


0.52 


53106 




A06 


SNLS-04D2mh 


09:59:45.872 


+02:08:27.94 


0.60 


53356 






SNLS-04D2mj 


10:00:36.535 


+02:34:37.44 


0.51 


53356 






SNLS-05D2ab 


10:01:50.833 


+02:06:23.02 


0.32 


53375 






SNLS-05D2ac 


09:58:59.244 


+02:29:22.22 


0.49 


53375 






SNLS-05D2bv 


10:02:17.008 


+02:14:26.05 


0.47 


53391 






SNLS-05D2cb 


09:59:24.592 


+02:19:41.34 


0.43 


53409 






SNLS-05D2dm 


10:02:07.611 


+02:03:17.35 


0.57 


53441 






SNLS-05D2dw 


09:58:32.058 


+02:01:56.36 


0.42 


53441 






SNLS-05D2dy 


10:00:58.083 


+02:10:59.52 


0.50 


53441 


2 




SNLS-05D2ec 


09:59:26.170 


+02:00:49.36 


0.53 


53441 


1,2 




SNLS-05D2ei 


10:01:39.103 


+01:49:12.02 


0.37 


53441 






SNLS-05D2hc 


10:00:04.574 


+01:53:09.94 


0.36 


53463 






SNLS-05D2ie 


10:01:02.907 


+02:39:28.90 


0.35 


53467 
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Table 2 — Continued 



Name 


"J2000.0 


<5j2000.0 


ZSPEC 


Discovery (MJD) 


Culling Status 3 - 


Refs b 


SNLS-04D3df 


14:18:10.020 


+52:16:40.13 


0.47 


53110 




A06 


SNLS-04D3ez 


14:19:07.916 


+53:04:18.88 


0.26 


53110 




A06 


SNLS-04D3fk 


14:18:26.212 


+52:31:42.74 


0.36 


53117 




A06 


SNLS-04D3gt 


14:22:32.594 


+52:38:49.52 


0.45 


53124 




A06 


SNLS-04D3hn 


14:22:06.878 


+52:13:43.46 


0.55 


53124 




A06,H05 


SNLS-04D3kr 


14:16:35.937 


+52:28:44.20 


0.34 


53147 




A06,H05 


SNLS-04D3nh 


14:22:26.729 


+52:20:00.92 


0.34 


53166 


1,2 


A06,H05 


SNLS-04D3nq 


14:20:19.193 


+53:09:15.90 


0.22 


53176 




A06,H05 


SNLS-05D3cf 


14:16:53.369 


+52:20:42.47 


0.42 


53410 






SNLS-05D3ci 


14:21:48.085 


+52:26:43.33 


0.51 


53416 






SNLS-05D3dd 


14:22:30.410 


+52:36:24.76 


0.48 


53441 






SNLS-05D3gp 


14:22:42.338 


+52:43:28.71 


0.58 


53462 






SNLS-05D3hq 


14:17:43.058 


+52:11:22.67 


0.34 


53474 






SNLS-05D3jq 


14:21:45.462 


+53:01:47.53 


0.58 


53474 






SNLS-05D3jr 


14:19:28.768 


+52:51:53.34 


0.37 


53474 


1.2 




SNLS-05D3kx 


14:21:50.020 


+53-08-1 3 49 


0.22 


53519 






SNLS-05D31q 


14:21:18.449 


+52:32:08.29 


0.42 


53528 






SNLS-05D3mq 


14:19:00.398 


+52:23:06.81 


0.24 


53559 


5 




SNLS-05D3mx 


14:22:09.078 


+52:13:09.35 


0.47 


53559 






SNLS-03D4ag 


22:14:45.790 


-17:44:23.00 


0.28 


52813 




A06 


SNLS-03D4au 


22:16:09.920 


-18:04:39.37 


0.47 


52815 




A06,B06 


SNLS-03D4cj 


22:16:06.660 


-17:42:16.72 


0.27 


52873 




A06,H05 


SNLS-03D4gf 


22:14:22.907 


-17:44:02.49 


0.58 


52930 




A06,B06 


SNLS-03D4gg 


22:16:40.185 


-18:09:51.82 


0.59 


52930 




A06,B06 


SNLS-03D4gl 


22:14:44.177 


-17:31:44.47 


0.57 


52935 


5 


A06,H05 


SNLS-04D4bq 


22:14:49.391 


-17:49:39.37 


0.55 


53174 




A06,B06 


SNLS-04D4gg 


22:16:09.268 


-17:17:39.98 


0.42 


53228 




H05 


SNLS-04D4gz 


22:16:59.018 


-17:37:19.02 


0.38 


53235 






SNLS-04D4ht 


22:14:33.289 


-17:21:31.33 


0.22 


53254 


2 




SNLS-04D4in 


22:15:08.585 


-17:15:39.85 


0.52 


53267 






SNLS-04D4jr 


22:14:14.335 


-17:21:00.93 


0.48 


53284 






SNLS-04D4ju 


22:17:02.733 


-17:19:58.34 


0.47 


53284 







a Numbers indicate the criteria from J3]that caused the SN to be rejected from the sample 

b References indicated as follows: A06 - lAstier et ah! j2006h . H05 - iHowell et"aH ^20051) . B06 - iBasa et all ll2006h 
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Table 3. SN la Sample 







Total 


Culled 


Field 


Season 






Dl 


1 


n 


8 




2 


9 


8 


D2 


1 


8 


4 




2 


13 


11 


D3 


1 


8 


7 




2 


11 


9 


D4 


1 


6 


5 




2 


7 


6 


ALL 




73 


58 



Table 4. Unconfirmed SN la Candidates 



Name QJ2000.0 <5j2000.0 zphot Discovery (MJD) Type XsNla Status 



SNLS-03Dlgc 


02: 


:24: 


06 


.043 


-04: 


:23: 


:19 


.14 





.54 


52993 


SN? 


4. 


.208 


Probable SN la 


SNLS-04D21u 


10: 


:01 


:09 


.465 


+02: 


32: 


14. 


52 





.37 


53353 


SN? 


3 


.741 


Probable SN la 


SNLS-04D21x 


10: 


01 


17 


.159 


+01: 


42: 


50. 


97 





.50 


53353 


SN? 


2 


.031 


Probable SN la 


SNLS-04D3M 


14: 


:16 


17 


.101 


+52: 


19: 


28. 


40 





.53 


53135 


SN? 


3 


.129 


Probable SN la 


SNLS-05D3ba 


14: 


:18 


:26. 


790 


+52: 


41: 


50 


56 





.44 


53387 


SN? 


2 


.323 


Probable SN la 


SNLS-05D31c 


14: 


:22: 


:22. 


.902 


+52: 


'28: 


:44. 


11 





.49 


53519 


SN 


2 


.027 


Probable SN la 


SNLS-05D31x 


14: 


:17 


56 


.809 


+52: 


20: 


23. 


26 





.58 


53532 


SN? 


3 


.102 


Probable SN la 


SNLS-03D4bx 


22: 


:14: 


:48. 


602 


-17: 


:31 


:17. 


.58 





.54 


52843 


SN? 


2. 


.089 


Probable SN la 


SNLS-03D4ev 


22 


:16 


:51 


.395 


-17: 


:20 


02 


.37 





.53 


52914 


SN? 


2. 


.591 


Probable SN la 


SNLS-04D4cm 


22 


:13: 


:28 


.782 


-18: 


:03 


40 


.56 





.55 


53177 


SN? 


3 


.928 


Probable SN la 


SNLS-04D4et 


22 


11 


:51 


788 


-17: 


:47: 


:22 


86 





.58 


53204 


SN? 


1 


.558 


Probable SN la 


SNLS-04D4iy 


22 


17 


:07. 


977 


-18: 


:07 


:07. 


18 





.51 


53267 


SN? 


2 


.801 


Probable SN la 


SNLS-04Dlfh 


02: 


:26 


:59. 


.401 


-04: 


29 


:42. 


.41 





.56 


53235 


SN? 


5 


.762 


Possible SN la 


SNLS-04Dlqn 


02: 


:27 


:28. 


186 


-04: 


20 


35 


78 





.52 


53323 


SN? 


6 


.400 


Possible SN la 


SNLS-04Dlqr 


02: 


:25: 


:49. 


.083 


-04: 


:29: 


:00 


.23 





.58 


53323 


SN? 


5 


.135 


Possible SN la 


SNLS-04Dlsc 


02: 


:26: 


:34. 


.371 


-04: 


:02 


:45 


60 





.59 


53351 


SN? 


6 


.038 


Possible SN la 


SNLS-05D3bz 


14: 


17 


:50 


119 


+52: 


.51: 


24. 


16 





.52 


53410 


SN? 


5. 


.499 


Possible SN la 


SNLS-05D3ch 


14: 


:19: 


:09 


.668 


+52: 


:47: 


35. 


93 





.35 


53416 


SN? 


5 


.553 


Possible SN la 


SNLS-04D4ay 


22: 


:15: 


:54. 


.038 


-18: 


:02: 


:48. 


.95 





.58 


53174 


SN? 


5 


.064 


Possible SN la 
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Table 5. Spectroscopic Completeness 



Field 




SNe la 




Completeness Fraction a 


(Inn firm pH 


Pro n^ nl p 


Po^i nip 


Minimum 


Most Likely (Cspec) 


Dl 


16 


1 


4 


0.76 


0.94 


D2 


15 


2 





0.88 


0.88 


D3 


16 


4 


2 


0.73 


0.80 


D4 


11 


5 


1 


0.65 


0.69 


ALL 


58 


12 


7 


0.75 


0.83 



a The maximum completeness is 1.00 



Table 6. i' Frame Limit Equation Parameters 



IQgood 


IQbad 


E re f 


Sref 


L .5 ot 


(arcsec) 


(arcsec) 


(s) 


(DN) 


(mag) 


0.69 


1.06 


3641 


29.18 


24.5 2.22 



Table 7. Monte Carlo Efficiencies 



On-Field Yearly 



Field %' Detection Spec Spec (e yr ) 



Dl 


0.948 


0.612 


0.299 


D2 


0.981 


0.528 


0.217 


D3 


0.971 


0.629 


0.313 


D4 


0.979 


0.654 


0.310 
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Table 8. SNLS Type la SN Volumetric Rates 





fRAW 


fobs 


n +z h 


e 


^G,0.2<z<0.6 


rv 


Field 


(yr- 1 ) 


(yr- 1 ) 


(y^ 1 ) 


(degrees 2 ) 


(xlO 4 Mpc 3 ) 


(xlO~ 4 yr" 1 Mpc~ 3 ) 


Dl 


26.7 ± 6.7 


28.4 ± 7.1 


41.7 ± 10.4 


1.024 


106.0 


0.39 ± 0.10 


D2 


34.6 ± 8.9 


39.4 ± 10.2 


57.7 ± 14.9 


1.026 


106.2 


0.54 ± 0.14 


D3 


25.5 ± 6.4 


31.9 ± 8.0 


46.8 ± 11.7 


1.029 


106.5 


0.44 ± 0.11 


D4 


17.7 ± 5.3 


25.7 ± 7.8 


37.7 ± 11.4 


1.027 


106.3 


0.35 ± 0.11 


AVG C 


24.1 ± 3.3 


30.3 ± 4.0 


44.4 ± 5.9 


1.026 


106.2 


0.42 ± 0.06 d 



a rates after correcting for spectroscopic incompleteness 
b rates after correcting for time dilation 
c Poisson error weighted averages 
Statistical error only 



Table 9. Summary of Uncertainties 



Source 57v a Srr h 



Poisson ±0.06 ±0.020 

Luminosity Estimate 



Total Systematic 



+0.033 
-0.023 



+0.03 +0.010 
-0.08 -0.031 



Spec. Completeness 

Host Extinction ±0.10 ±0.037 

Frame Limits -0.03 -0.011 

Stretch ±0.01 ±0.004 

Subluminous SNe la ±0.08 ±0.029 



Total Statistical ±0.06 



+0.13 +0.048 
-0.09 -0.033 



a Volumetric uncertainty in units of 
10~ 4 yr- 1 Mpc~ 3 

b Luminosity specific uncertainty in 
SNu 



